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The  motivation  for  determining  the  response  time  for  a  pressure  regulator 
to  adjust  the  pressure  at  an  area  of  interest  located  at  the  end  of  a  long  tube 
may  be  due  to  several  reasons.  One  of  these  is  the  possibility  that  the  regulator 
cannot  be  conveniently  controlled  from  a  more  desirable  position  closer  to  the 
area  of  interest.  Another  reason  to  place  the  regulator  a  long  distance  from  the 
area  of  interest  might  be  to  avoid  high  pressure  lines  over  a  prolonged  distance. 
Consider  the  system  shown  in  Figure  1,  which  is  the  proposed  propulsion 
system  for  a  microgravity  experiment  located  in  the  Shuttle  cargo  bay.  The 
regulator  is  controlled  electronically  so  the  ability  to  reduce  the  amount  of 
electrical  wiring  will  increase  system  reliability.  Furthermore,  Shuttle  safety 
requirements  demand  that  high  pressure  lines  be  tested  for  space  worthiness, 
thus  placement  of  the  regulator  close  to  the  propellant  source  may  eliminate  the 
presence  of  these  high  pressure  lines. 

10  to  0.3  MPa  variable  regulator 


Figure  1:  Proposed  Microgravity  Propulsion  System 
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To  determine  the  time  delay  for  the  system  above,  it  is  assumed  that  the 
pressure  loss  due  to  gas  escaping  out  the  nozzle  throat  area  is  negligible  and 
that  the  initial  start-up  delay  will  be  the  longest  delay  encountered.  The  latter 
assumption  is  valid  when  one  considers  that  the  pressure  in  the  tube  and 
plenum  are  essentially  zero  before  the  regulator  is  activated.  For  a  worst  case 
scenario,  we  will  assume  that  the  first  command  to  the  regulator  will  demand  the 
maximum  pressure  of  10  MPa. 

Given  a  demand  of  10  MPa  in  the  plenum,  the  first  step  in  this 
pneumatics  problem  is  to  determine  the  effective  pipe  length  and  pressure  loss 
in  the  system.  Assuming  30  ft  of  ABS  (thermoplastic  polymer)  or  PVC  pipe  is 
used  and  that  four  90°  elbow  joints  are  needed,  the  effective  pipe  length  may 
be  found  from: 

Leff  *=  L$t+  nCFD1*2 

where  L  is  the  length  of  straight  pipe,  n  is  the  number  of  joints,  C  is  a  constant 

equal  to  43.7  for  length  in  feet  and  diameter  in  inches,  F  is  the  resistance 
coefficient  of  the  fittings  used  (0.5  for  90°  elbow  joints),  and  D  is  the  pipe 
diameter.1  For  the  system  shown,  the  pipe  diameter  will  be  1/4  inch  or  less. 

The  pressure  loss  associated  with  a  given  effective  length  is  found  from: 

* 

AP  =  FLQ  /  PD*-33 

where  AP  is  the  pressure  loss  in  bars,  F  is  a  friction  factor  (500  for  steel  pipe, 
variable  for  ABS  or  PVC  pipe),  L  is  the  effective  length,  Q  is  the  flow  rate,  P  is 
the  final  pressure  (10  MPa,  or  100  bar),  and  D  is  the  pipe  diameter  in  mm.1 
These  figures  result  in  a  calculated  pressure  loss  of  about  3  bar  (0.3  MPa),  or 
only  3%  of  the  desired  pressure,  for  a  pipe  diameter  of  1/4  in  (6.35  mm)  and  a 
flow  of  27.43  liters/sec.  As  a  result,  it  is  safe  to  assume  negligible  pressure  loss. 


Having  determined  the  effective  pipe  length  and  pressure  loss,  the  delay 
time  for  the  pipe  and  plenum  to  pressurize  to  the  desired  10  MPa  may  be 
determined.  Assuming  that  the  gaseous  nitrogen  acts  as  a  perfect  gas,  and  that 
conservation  of  mass  momentum,  and  energy  hold,  we  may  use  the  following 
set  of  coupled  equations:2 


Perfect  Gas  Law: 

P  =  pRT 

Momentum  Equation: 

Tdu  dul  dP  _ 

K  [dt  dHj  d« 

Continuity  Equation: 

B,K,^^,pBu)'° 

Energy  Equation: 

P -Me.  v)*— r 

dt '  2  9  da  L 

p  (E  +  Puj  -  0 

where  P  is  pressure,  p  is  density,  R  is  the  gas  constant  for  GN2  (296.7972 

m2/s2K°),  T  is  absolute  temperature,  v  is  the  gas  velocity,  A  is  the  pipe  cross 
sectional  area,  E  is  the  product  of  the  specific  volume  and  temperature,  and  x  is 

4 

the  location  down  the  pipe  (x  =  0  at  the  regulator,  x  =  L  at  the  plenum). 

Assuming  the  temperature  of  the  gas  remains  constant  at  -20°  C  (253.15  °K), 
the  energy  term  (CvT)  remains  constant.  As  a  result,  the  energy  derivatives  are 

zero  and  the  associated  equation  reduces  to  the  momentum  equation  leaving 
three  equations  and  three  unknowns  (P,  p,  v).  These  equations,  however,  are 
coupled,  presenting  a  nontrivial  problem  when  solved  analytically.  For 
simplicity,  then,  a  numerical  solution  to  this  set  of  equations  will  be  sought, 
assuming  the  errors  involved  are  acceptable. 


Given  that  the  regulator  valve  opening  time  may  be  considered 
instantaneous  (typical  valve  opening  times  are  around  0.01  to  0.03  sec.),  we 
may  assume  that  the  effective  pipe  cross  sectional  area  is  a  constant,  so  the 
force  acting  on  the  gas  is  also  a  constant  (F  =  PA).  Given  this  assumption,  the 
time  for  Poiseuille  flow  to  form  within  the  pipe  may  be  found.  This  flow 
represents  the  parabolic  flow  of  a  fluid  or  gas  in  a  pipe  due  to  a  pressure 
gradient.  Once  Poiseuille  flow  has  formed,  the  fluid  traveling  along  the 
centerline  of  the  pipe  is  at  a  maximum  velocity  for  the  given  pressure 
differential.  For  starting  flow  in  a  pipe,  it  may  be  assumed  that  the  time  for 
Poiseuille  flow  to  form  is  given  by:2 

t  =  0.75r2p/p 

where  r  is  the  pipe  radius,  p  is  the  gas  density,  and  p  is  the  gas  viscosity.  Given 
that  GN2  at  10  MPa  and  -20°C  has  a  density  of  154  kg/m3  and  a  viscosity  of 

14.826  kg/m  sec'1,  the  Poiseuille  flow  in  a  1/4  inch  diameter  pipe  occurs  in  8  x 
10'5  sec,  which  may  be  assumed  to  be  instantaneous. 

Given  that  the  development  of  Poiseuille  flow  will  not  cause  a  significant 
time  delay,  the  time  for  the  pipe  and  plenum  to  pressurize  may  be  found  by 
determining  the  mass  flow  rate  of  the  gas  and  the  associated  rate  of  pressure 
chdnge.  Next,  performing  numerical  integration  by  taking  infinitesimally  small 
time  steps,  the  pressure  in  the  pipe  and  plenum  may  be  found.  Finally,  this 
pressure  is  used  to  determine  a  new  mass  flow  rate  and  the  process  is 
continued  until  the  desired  pressure  is  achieved. 

Four  assumptions  associated  with  this  process  are:  1)  the  pipe  has  some 
very  small  initial  pressure,  2)  the  change  in  pressure  within  the  pipe  and 
plenum  is  uniform  down  the  entire  pipe  and  plenum  length,  3)  the  friction  loss  in 
the  pipe  is  negligible  due  to  the  earlier  calculation  of  the  effective  length,  and  4) 


the  change  in  fluid  velocity  over  each  time  step  (dv/dt)  is  zero. 

Determining  the  mass  flow  rate  involves  dividing  the  momentum 

equation  by  pv2  and  replacing  p  with  P/RT  and  v2  with  (m/Ap)2  to  yield: 


dp  +  du  =  0 
mzRT  v 


integrating  each  term  yields: 


reg 


solving  for  the  mass  flow  rate: 


where  Pr  is  the  pressure  leaving  the  regulator  (10  MPa),  P  is  the  pipe  and 

reg 

plenum  pressure  and  is  initially  not  zero,  and  A,  R,  and  T  are  as  previously 
defined. 

Differentiating  the  perfect  gas  law  and  replacing  the  density  with  the  fluid 
mass/pipe  and  plenum  volume  (V),  produces: 

P  =  mRT/U 


Having  determined  the  rate  of  pressure  change  for  the  given  time  step, 
the  pipe  and  plenum  pressure  may  now  be  found  from  a  simple  Euler  equation: 


new 


»  P 


old 


Pdt 


Finally,  updating  the  time  so  t  =  t  +  dt  completes  the  cycle.  Repeating  this 
pattern  over  infinitesimally  small  time  steps  until  P  =  Pf0g  will  yield  the  time  to 


pressurize  the  pipe  and  plenum.  Using  this  approach,  the  time  to  pressurize  the 
system  from  0  to  10  MPa  for  various  pipe  diameters  is  shown  in  Figure  2. 

Once  the  system  has  been  pressurized  to  10  MPa,  the  next  concern  may 
be  the  time  delay  involved  in  reducing  the  system  pressure,  as  opposed  to 
building  it  up.  For  the  particular  system  shown  in  Figure  1 ,  the  only  way  to 
reduce  the  pressure  is  by  lowering  the  pressure  output  from  the  regulator  and 
letting  gas  bleed  from  the  nozzle.  A  modification  to  the  system  might  include  a 
relief  T-valve  on  the  system  to  vent  the  excess  gas.  Although  this  would  cause 
a  discontinuity  in  microgravity  acceleration  measurements,  the  pressure 
reduction  would  be  nearly  instantaneous.  To  determine  the  bleed-down  time 
associated  with  just  the  nozzle,  the  same  numerical  integration  approach  used 
earlier  applies.  One  major  difference,  however,  is  the  mass  flow  rate  equation 
is  for  a  supersonic  nozzle  and  is  a  function  of  the  throat  area:3 


m 


fltPk 


(2/dc+D) 


(k+n/oc-n 


kRT 


where  k  is  the  ratio  of  specific  heats  which,  for  nitrogen  operating  in  a  pressure 

range  of  0.3  to  10  MPa  at  -20°C,  is  1 .432  to  1 .783.  Also,  the  updated  pressure 
•  • 

is  Pnew  =  poid  *  Pdt-  Fi9ure  3  shows  the  time ‘delay  for  the  system  to  bleed  from 
10  to  0.3  MPa  given  various  pipe  diameters. 

.  In  conclusion,  a  numerical  solution  to  determine  the  time  delay 
associated  with  pressurizing  an  area  of  interest  due  to  placing  the  pressure 
regulator  upstream  from  that  area  has  been  presented.  Although  several 
assumptions  have  been  made,  by  taking  a  worst  case  scenario  approach  the 
resulting  solution  should  yield  a  rough  idea  of  the  anticipated  time  delay. 


pressure  (MPa)  pressure  (MPa) 


Pressure  Build-Up  Times  for  Various  Pipe  Diameters 


Pressure  Bleed  Times  for  Various  Pipe  Diameters 


time  (sec) 

Figure  3 


References 


1.  Waning,  R.H.  Pneumatic  Handbook.  6th  ed.  Gulf  Publishing, 
Houston,  1982.  Pp  146-155. 

2.  White,  Frank  M.  Viscous  Fluid  Flow.  Mcgraw-Hill,  1974. 

Pp  112,  140-145,  661. 

3.  Sutton,  George  P.  Rocket  Propulsion  Elements.  5th  ed. 

Wiley  &  Sons,  New  York,  1986.  Pp  46. 


Determination  of  Shuttle  Limit  Cycles  for 
Various  Angles  of  Attack  and  Deadbands 

by  Scott  Baune 


The  purpose  of  this  paper  is  to  analytically  determine  the  limit  cycle  of  the 
shuttle  about  the  pitch  axis  for  various  angles  of  attack  and  deadbands.  The 
angle  of  attack,  a,  is  defined  by  the  rotation  about  the  shuttle  y-axis  so  that  a  = 
90°  is  gravity  gradient  stable.  The  deadband  is  the  region  of  deviation  from  a 
that  the  shuttle  will  be  allowed  to  drift  before  corrections  are  implemented  by  the 
Reaction  Control  System  (RCS).  Therefore,  an  a  of  50°  with  a  deadband  of  ± 

5°  would  allow  the  shuttle  to  drift  from  a(  =  45°  up  to  ah  =  55°  without  RCS 
control. 

The  primary  source  of  disturbance  about  the  pitch  axis  above  an  altitude 
of  200  km  is  assumed  to  be  the  gravity  gradient  torque  given  by: 

^moK*^  2a 

where  is  defined  as  the  maximum  gravity  gradient  torque  for  a  given 
altitude.  This  torque  occurs  at  a=rc/4  and  is  given  by: 


P 

max 


3GMe 
2(h  +  Re)  3 


where  GMa  =  398,601 .2  km3/s2,  Re  =  6378.137  km,  lxx  =  1 ,340,000  Nms2,  lzz  = 
10,200,000  Nms2,  and  h  is  the  orbital  altitude. 

To  determine  the  limit  cycle  due  to  gravity  gradient,  the  moments  acting 
about  the  pitch  axis  are  summed  as: 


IM  -  I y y  (d2a/dt2)  -  P,„0II  «ln2a 

where  lyy  =  9,910,000  Nms2.  It  should  be  noted  that  the  given  moment  of 
inertia  values  are  typical  but  will  vary  from  mission  to  mission. 

At  this  point  the  notation  shown  in  Figure  One  is  implemented  so  that  0  + 
a  =  900  and  (d2o/dt2)  =  -(d2a/dt2). 


0° 

90 


a  -  90  0 

<f>=  0° 


Figure  One:  Ongle  Notation 


As  a  result  of  this  change,  the  moment  equation  may  now  be  written  as: 

2M  =  lyy  (d2er/dt2)  -  -Pm8„  *in2* 


Integration  of  this  equation  produces: 


(dar/dt)2  -  CPmoii/lyy)  co$2a  *  C 


where  C  is  a  constant  of  integration  found  from  the  initial  conditions  that  de/dt 
and  0  at  t  =  0  are  0  rad/sec  and  eh>  respectively.  Therefore,  substitution  of  C 
and  implementation  of  cos2e  =  1  -  2sin2o  yields: 

(df/dt)2  -  2(PmaM/lyy)  *ln2*h  (1  -  n2sln2») 

where  n  =  1/sin0h. 

Solving  for  the  time  to  drift  from  some  angle  a  to  gravity  gradient  stable  (a 
=  90°  and  a  =  0°)  yields: 


t  = 


N 


2  man  sin  2<|>t 
1  aa 


I 


d<j> 


NJ 


1  -  n  sin  2$ 


By  definition,  the  elliptic  integral  above  may  be  stated  as: 

t  -  kc  F(fl,k) 

where  k  is  known  as  the  modulus  and  is  given  by  1/n,  or  sinoj,,  c  is  the  constant 
term  in  front  of  the  integral  above,  and  F(B,k)  is  the  value  of  the  elliptic  integral 
easily  found  in  tabulated  form  or  computationally  with  B  given  by: 


sinfih>|  -  n  *ln»h  , 


To  determine  the  total  drift  time,  T,  across  a  region  from  0|,  through  e*,, 
and  back  to  e>\  as  shown  in  Figure  Two,  the  time  difference  between  these  two 
integrals  must  be  found,  then  doubled.  The  time  is  doubled  because  the 
difference  between  the  two  integrals  from  a ^  to  Q\  only  accounts  for  the  upper 
half  of  the  curve  in  Figure  Two  based  on  the  given  initial  conditions.  Thus  the 
total  drift  time  is  given  by: 


T  -  2kc  { F ( B h , k )  -  F(fl„k)} 


As  shown,  Figure  Two  illustrates  the  the  drift  rate  of  the  shuttle  about  the 
pitch  axis  versus  the  pitch  angle.  The  initial  conditions,  that  is,  do/dt  =  0  rad/sec 
and  0  =  0f,  at  t  =  0,  used  in  the  elliptic  integral  formulation  may  then  be  seen  at 
the  left  portion  of  the  curve,  as  well  as  the  pitch  rate  at  the  end  condition,  a\ . 
Finally,  returning  to  the  the  lower  position  on  Figue  Two  requires  the  use  of 
RCS  thrusters  and  is  considered  instantaneous  in  the  formulation  above. 
However,  if  the  force  of  the  thrust,  F,  and  the  distance  from  the  shuttle  center  of 
gravity,  0,  are  known,  the  thrust  duration  may  be  found  from: 

t  -  (2lyy/FD)  der/dt 

Implementation  of  this  derivation  was  accomplished  usinfi  an  altitude  of 
275  km,  where  Pmax  =  18  nm.  Figure  Three  shows  the  drift  time  in  minutes  for 
various  angles  of  attack  and  deadbands  of  ±  2,  3,  4,  5,  and  10°. 

To  check  the  results  of  this  graph,  a  new  formulation  which  holds  the 
gravity  gradient  torque  constant  at  its  maximum  level,  despite  changes  in  a,  was 
made  so  that: 


»(t)  -  -(  (PmaH/2lyy>  **n2ar  >  t2  +  *h 


and  the  total  drift  time,  T,  is  given  by: 


T  =  2t 


M 


4>h  -  4>i 


Ijujul  sin  2<t> 

21  yy 


Figure  Four  compares  the  results  of  the  fixed  torque  formulation  to  Figure 
Three  for  deadbands  of  ±  2  and  10°.  This  comparison  reveals  similiar  results 


around  the  a  =  45°  area  but  not  around  the  ends.  These  results  appear 
reasonable  since  the  gravity  gradient  torque  at  45°  is  the  same  for  the  fixed 
torque  case  and  the  elliptic  integral  case,  but  as  a  approaches  0  or  90°,  the 
gravity  gradient  torque  reduces  to  zero  in  the  elliptic  integral  case,  while 
remaining  constant  in  the  fixed  torque  case. 


aoa  (deg) 


Figure  Three 


Drift  Time  vs.  AOA  for 

+/-  10  degree  deadband 


Drift  Time  vs  AOA  for 

+/-  2  degree  deadband 
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SCOTT  A.  BAUNE 


Design  of  a  Variable  Thrust  Cold  Gas  Microthruster  Nozzle 

The  motivation  for  the  design  of  a  variable  thrust  microthruster  stems  from 
the  desire  to  counter  the  gravity  gradient  and  drag  accelerations  present  on  a 
microgravity  facility  or  experiment  located  at  positions  other  than  the  center  of 
gravity  of  another  facility  (i.e.  the  Shuttle  or  Space  Station).  Research  in  the 
development  of  microthrusters  began  as  early  as  the  1950’s  with  the 
development  of  hot  gas  thrusters.  These  thrusters  were  scaled  versions  of 
larger  variable  thrust  engines  which  had  already  been  developed.  Failure  of 
these  thrusters  was  mainly  a  result  of  materials  failure  caused  by  the 
combustion  of  the  propellant.  Recently,  however,  microthrusters  using  a  cold 
gas  propellant  such  as  nitrogen  or  ammonia  have  been  developed  for  attitude 
control  and  microgravity  research. 

These  systems  incorporate  a  regulator  to  vary  the  pressure  of  the  gas 
inside  a  collection  chamber.  The  gas  then  passes  through  a  nozzle  via  a 
constricted  region  known  as  the  throat,  producing  thrust  at  the  nozzle  exit.  The 

design  of  the  nozzle  for  a  cold  gas  consists  of  determining  certain  critical 
parameters  shown  in  Figure  1.  These  include  throat  area,  Af,  chamber 

pressure,  P1f  chamber  temperature,  T  ,  divergence  angle,  a,  and  expansion 
ratio,  Ag/Aj. 

Assuming  that  isentropic  relations  hold  inside  the  chamber  and  nozzle 
generally  allows  us  to  calculate  performance  to  within  10%  of  the  measured 
values.  As  a  result,  the  design  process  illustrated  here  is  limited  to  this  ideal 
rocket  case.  Once  the  ideal  rocket  analysis  is  complete,  the  results  must  be 
corrected  to  account  for  the  following  actual  conditions  present  within  the 
engine:1 


1 .  Conical  divergence  of  the  exhaust  jet. 

2.  Surface  friction  and  flow  disturbances  in  the  exhaust  nozzle. 

3.  Constrictic  i  of  the  exit  area  due  to  boundary  layer  build-up. 

4.  Jet  detachment. 

5.  The  gas  velocity,  pressure,  temperature,  and  density  are  not  uniform 
across  any  section  normal  to  the  nozzle  axis. 

6.  The  propellant  is  not  a  true  perfect  gas  and  is  not  homogeneous  and 
invariant  within  the  chamber  and  nozzle.2 

The  last  condition  allows  the  implementation  of  the  perfect  gas  law. 

0)  v  p»rVm 


where  p,  p,  and  T  are  the  pressure,  density,  and  temperature,  respectively,  at 
some  point  x,  R  is  the  universal  gas  constant  (8314.3  Joule/k  mol-°K),  and  M  is 
the  molecular  mass  of  the  propellant.  For  a  perfect  gas,  the  specific  heat  at 


constant  pressure,  cp,  the  specific  heat  at  constant  volume,  cyl  and  the  ratio  of 

the  two  are  constant: 

(2)  k  -  c  /c 

p  u 

cp  "  cw  “  R/(MJ> 

(3»  CP- 

where  J  =e  1  in  SI  units  and  778  ft-lb/Btu  in  English  units. 

Applying  the  principle  of  conservation  of  energy  in  which  the  decrease  in 
enthalpy  is  equal  to  the  increase  in  kinetic  energy  of  the  propellant,  we  may 
write  the  enthalpy,  which  is  the  product  of  the  specific  heat  and  the  absolute 
temperature,  as: 

(4)  ht  -  h2  -  cpttj  -  T2)  -  <u22  -  u, 2)/2J 

where  h  is  the  enthalpy  and  v  is  the  velocity  at  the  chamber  and  the  nozzle  exit. 
Solving  (4)  for  the  nozzle  exit  velocity,  v2,  yields: 

(5)  uz  ■  VV*".2 

In  addition,  for  isentropic  flow,  the  following  relations  hold2: 

(6)  T,/T2  >  (p|/p2)(,‘-',/lt  -  (p2/P|) IM 


When  the  chamber  cross  section  is  large  compared  to  the  nozzle  cross  section, 
the  chamber  velocity,  vr  is  small  and  may  be  neglected.  Thus  (5),  (6),  and  (3) 

may  be  written  as: 


(7) 


Given  the  exit  velocity  it  is  now  possible  to  determine  the  specific 
impulse,  I  ,  of  the  fuel: 

*r 

-  ^<n>g0) 

where  F  is  thrust,  m  is  the  mass  flow  rate  of  the  propellant,  and  gQ  is  standard 
sea-level  gravity  (9.80  m/s2).  The  thrust  may  be  defined  as  the  integral  of  the 

pressure  forces  acting  on  the  nozzle  and  chamber  in  the  axial  direction.  This 

integral  equates  to: 

F  -  mu2+  <p2-  p3)f»2 

where  the  first  term  is  known  as  the  momentum  thrust  and  the  second  as  the 
pressure  thrust.  However,  if  the  exit  pressure  is  assumed  to  be  equal  to  the 
ambient  pressure,  the  pressure  thrust  term  is  zero  and  the  specific  impulse  may 
be  written  as: 

<8>  *  U2/9o 

To  determine  the  thrust,  the  mass  flow  rate  is  needed.  The  principle  of 
conservation  of  matter  is  used  and  is  expressed  by  equating  the  mass  flow  rates 
at  any  two  sections: 

•  « 

(9)  mi“m2-fl1u1pJ-fl2u2p2 

Applying  (6)  and  (7)  to  (9)  yields  the  propellant  mass  flow  rate  in  any  section  as 
a  function  of  the  chamber  temperature  and  pressure  ratio: 


(10) 


The  maximum  propellant  flow  per  unit  area  occurs  at  the  throat,  so  the  throat 
pressure,  p ,  may  be  found  by  differentiating  (10)  with  respect  *o  x  and  setting 

the  derivative  to  zero: 


(11) 


2 

k+1 


k/0e-n 


From  (11)  it  is  now  possible  to  determine  the  conditions  for  which 
subsonic,  sonic,  and  supersonic  flow  at  the  throat  occur.  Referring  to  Figure  1 , 

sonic  or  supersonic  conditions  at  the  throat  occur  if  the  chamber  to  exit  pressure 
ratio  (pt/p2)  is  equal  to  or  greater  than  p(/pr  Subsonic  throat  conditions  occur  if 

p1/p2  is  less  than  p/pr  In  addition,  chamber  pressures  below  a  given  level  will 


also  result  in  subsonic  flow.  This  level  is  approximately  2.17  atm  for  nitrogen 
and  2.6  for  ammonia  when  operating  at  sea-level.  Typically,  the  supersonic 
nozzle  is  used  for  rocket  engines  because  the  lower  velocities  associated  with 
subsonic  flow  tend  to  lower  the  specific  impulse  below  practical  levels  for  space 
flight. 


Continuing  our  determination  of  the  rocket  thrust,  we  combine  (11)  with 
the  isentropic  relationships  in  (6)  to  determine  the  density  and  temperature  at 
the  throat: 


(12) 


2 

k+1 


t/(k-n 


and 


(13) 


T  -  ill 

*  k+1 

and  combining  (7),  (11),  and  (13)  results  in  the  throat  velocity  equation: 


(14) 


u. 


2kRT, 
M(k+1 ) 


Thus  we  may  now  write  the  mass  flow  rate  for  a  supersonic  nozzle  as  a  function 
of  those  variables  which  are  design  controlled:  throat  area,  chamber  pressure, 
chamber  temperature,  and  specific  heat  ratio: 


(15) 


m  »  Atp,k 


(2/<k+U)0M)/0c'u 

kRT, 


Likewise,  using  the  perfect  gas  law,  and  equations  (7),  (12),  and  (14),  it  is 

possible  to  relate  another  important  design  parameter,  the  expansion  or  area 
ratio,  A^,  to  the  pressure  ratio  PjPy 


(16) 


' 

le+ 1 

l/(k-l) 

P 

t/k\  /  ic+t  r  p. 

k-l/k" 

2 

Pi 

v  J 

V  k’1  pi 

La 

Also,  given  the  expansion  ratio,  the  throat  area,  and  the  divergence 
angle  of  the  nozzle,  a,  the  nozzle  length  may  be  derived  from  simple  geometry: 


(17) 


rt  tan2(«) 


Finally,  we  are  in  a  position  to  relate  the  thrust  ot  a  supersonic  nozzle  to  the 
design  parameters.  Given  that  the  thrust,  defined  earlier,  may  be  wntten  as: 

F  -  mu2  +  (p2  -  p3)fl2 

Note  that  in  orbit  the  external  pressure  p3  is  negligible,  the  thrust  equation 
above  simplifies  to: 

F  =  mu2  +  p2fl2 

It  is  also  interesting  to  note  that  the  maximum  thrust  occurs  when  p2  is 

zero,  causing  V2  to  be  maximized.  Once  again  by  combining  equations  (7), 
(12),  and  (14)  we  may  relate  the  thrust  to  the  design  parameters  A(,  pr  k,  and 

A2/A(  (which  gives  Pg/p^  in  a  form  known  as  the  ideal  thrust  equation: 


(18) 


F 


l<k*t) 
_1_  <k-i) 

k+1 


(k-D/k 


*  fl2p 


2 


Another  form  of  the  ideal  thrust  equation  normalizes  the  force  by  dividing 
(18)  by  the  throat  area  and  chamber  pressure,  thus  the  unitless  variable  F/AtPi 
becomes  a  function  of  the  area  ratio  A2/A(  and  the  pressure  ratio  P2/P1  only. 


Validation  Case 


To  evaluate  the  above  equations,  data  from  a  cold  gas  ammonia 
propulsion  system  designed  for  the  Lincoln  Experimental  Satellites  LES-8  and 
LES-9  were  used.3  The  propulsion  system  for  both  satellites  was  provided  by 
TRW  Systems  Group  under  contract  to  M.l.T.  Lincoln  Laboratory.  Both  have 
identical,  cold  ammonia  propulsion  systems  with  a  specific  impulse  of  100 
seconds  operating  in  a  design  range  of  -10°C  to  40°C.  Given  the  15° 
divergence  angle  and  an  expansion  ratio  of  60:1,  a  nozzle  length  of  0.746  cm 
may  be  calculated.  Finally,  with  a  1.55  x  10'7  m2  throat  area  and  a  chamber 
pressure  of  about  0.207  MPa,  thrust  levels  of  0.045  N  were  achieved. 

Applying  this  data  to  the  equation  presented  earlier,  it  is  possible  to 
construct  a  plot  of  I  versus  area  ratio  as  shown  in  Figure  2  and  to  confirm  an 

I  of  100  seconds  for  an  area  ratio  of  60:1.  Likewise,  Figure  3  illustrates  a  plot 

of  equation  (16)  from  which  one  may  determine  that  the  pressure  ratios 
produced  from  a  60:1  area  ratio  at  chamber  temperatures  of  -10°C  and  40°C 
are  about  950  and  1725,  respectively.  These  ratios  represent  the  minimum  and 
maximum  pressure  ratios  obtainable  for  an  ammonia  system  having  the  given 
expansion  ratio  operating  over  the  given  temperature  range.  Finally,  the  thrust 
is  obtained  from  equation  (18).  Here  we  see  the  theoretical  thrust  levels 
obtainable  from  a  supersonic  nozzle  with  a  throat  area  of  1.55  x  10  m2,  a  60:1 
expansion  ratio,  a  chamber  pressure  of  0.207  MPa,  and  pressure  ratios  of  950 
to  1725  are  0.056  to  0.060  Newtons. 

The  disparity  between  the  calculated  thrusts  and  the  Lincoln  measured 
thrust  of  0.045  Newtons  is  a  result  of  the  ideal  nozzle  assumptions  made  earlier. 
The  first  of  these  accounts  for  conical  divergence  of  the  exhaust  jet.  in 
comparison  to  the  ideal  nozzle  flow,  which  assumes  axially  aligned  flow,  the 


Isp  vs.  Area  Ratio  for  Ammonia 

(PI  =  207000  Pa,  T1  =  -10  to  40  deg  C) 


Area  Ratio  Vs.  Pressure  Ratio 


for  ammonia  (T1  =  - 10  lo  40  deg  C) 


Figure  3 


real  rocket  contains  some  element  of  non  axially  directed  flow  which  results  in  a 
thrust  loss.  This  loss,  given  in  Altman,  may  be  simplified  to: 

F  -  Ffl(1  *  cos  a)  /  2 

where  Fq  is  the  ideal  thrust,  and  a  is  the  divergence  angle.  Thus  for  an  a  of  20°, 

the  thrust  due  to  conical  divergence  is  about  97%  of  the  ideal  thrust.  Table  1 
lists  the  thrust  correction  factor  due  to  various  divergence  angles. 


Divergence 

Thrust 

Angle 

Correction 

(deg) 

Factor 

0 

1.0000 

10 

0.9924 

20 

0.9698 

30 

0.9330 

Table  1:  Thrust  Correction  Factors 

for  ilarious  Divergence  Angies2 

The  second  ideal  rocket  assumption  neglected  surface  friction  and  flow 
disturbances.  Losses  due  to  these  effects  are  difficult  to  calculate,  however, 
they  may  be  minimized  by  reducing  sharp  curves  at  the  throat.  Typical  losses 
due  to  friction  flow  are  about  1  to  3%. 

The  third  assumption  neglected  constriction  of  the  exit  area  due  to 
boundary  layer  effects.  These  effects  are  most  prominent  at  the  nozzle  exit  and 
have  the  effect  of  reducing  the  operating  expansion  ratio.  According  to  Altman, 
observations  of  boundary  layer  growth  in  the  Reynolds  number  region  from  5  x 
105  to  5  x  10®  and  a  Mach  number  range  of  1.5  to  2.5  produce  the  following 
approximate  relationship  between  boundary  layer  displacement  thickness,  d, 


and  nozzle  length,  L: 


d  =  0.004L 

Applying  this  equation  and  (17)  to  the  Lincoln  example,  the  nozzle  length  is 
0.0041  m,  assuming  a  divergence  angle  of  20°,  and  the  boundary  layer 
displacement  is  1.65  x  10  m.  Simple  geometry  reveals  that  the  expansion 
ratio  has  been  reduced  from  60:1  to  59.856:1.  However,  this  reduction  in 
expansion  ratio  only  produces  a  thrust  loss  of  0.005%. 

Likewise,  the  remaining  assumptions  have  negligible  effects  on  the  thrust 
output  and  are  difficult  to  calculate  outside  of  experimentation.  According  to 
Altman,  the  combined  thrust  loss  due  to  these  assumptions  is  less  than  5%. 

Having  applied  these  corrections  to  the  calculated  thrust,  a  15  to  20% 
difference  still  exists  between  the  calculated  and  measured  thrust  reported  for 
the  LES  system.  This  difference  is  unaccounted  for  but  may  be  due  to  errors  in 
the  measured  thrust  data,  calculated  data,  or  the  interpretation  of  the  LES 
system  design. 


Design  Process 


For  the  Shuttle,  the  accelerations  due  to  gravity  and  drag  are  the  primary 
disturbances  in  the  low  frequency  (<  0. 1  Hz)  range.  These  forces  may  vary  from 
0  to  5  Newtons,  depending  on  the  Shuttle  angle  of  attack.  In  order  to  counter 
these  forces,  a  0  to  2  Newton  microthruster  using  gaseous  Nitrogen  as  the  cold 
gas  propellant  has  been  designed.  Four  of  these  thrusters  strategically  placed 
within  the  cargo  bay  may  then  be  used  to  create  a  microgravity  environment  of 
10‘6  G  or  smaller  for  a  position  off-set  from  the  Shuttle  center  of  gravity. 

Nitrogen  was  the  gas  of  choice  due  to  its  low  molecular  weight  and 
density.  This  allows  for  high  specific  impulses  and  the  ability  to  store  more 
propellant.  Originally  a  two-phase  system  was  considered  because  of  its 
increased  propellant  storage  capability,  however  system  complexity  and 
difficulties  involved  in  phase  separation  soon  outweighed  this  benefit. 


Isp  Vs.  Area  Ratio  for  Nitrogen 


Using  gaseous  nitrogen,  and  assuming  an  average  cargo  bay 
temperature  (chamber  temperature)  of  -20°C  and  a  chamber  pressure  range 
of  0.3  to  10  MPa,  a  range  of  k  values  from  1 .432  to  1 .783  may  be  obtained  from 
chemical  properties  tables.  The  pressures  used  roughly  correspond  to  the 
minimum  pressure  allowable  for  supersonic  flow,  and  25%  of  the  maximum 
pressure  obtainable  from  a  flight  qualified  cold  gas  storage  tank  (3000-5000 
psi).  The  25%  was  chosen  as  being  the  lowest  storage  pressure  allowable  to 
maintain  the  desired  thrust  levels. 

Assuming  that  throat  areas  below  1  x  10*7  m2  require  the  use  of  fuel 
filters,  and  that  excessive  nozzle  size  and  mass  will  not  be  inhibitory,  we  are  left 
only  to  decide  the  expansion  ratio,  throat  area,  and  length  of  that  will  produce  0 
to  2  Newtons  of  thrust  for  the  given  pressure  range. 

In  order  to  maximize  the  specific  impulse  of  nitrogen,  the  lsp  versus  area 
ratio  has  been  plotted  in  Figure  4.  From  this  figure  we  see  that  an  expansion 
ratio  of  about  8:1  will  yield  95%  of  the  theoretical  maximum  lsp  for  a  k  value  of 
1.783,  while  a  ratio  of  about  64:1  will  yield  95%  of  the  maximum  lsp  for  a  k  of 
1.432.  As  a  result,  an  expansion  ratio  of  64:1  will  maximize  lsp  efficiency  for 
either  value  of  k. 

Given  an  expansion  ratio  of  64:1,  equation  (16)  may  now  be  used  to 
determine  the  resulting  pressure  ratio  (P1/P2)  associated  with  that  area  ratio. 
Figure  5  shows  the  relation  between  the  two  ratios  and  reveals  that  for  a  64:1 
expansion  ratio,  a  pressure  ratio  of  2355  will  be  achieved. 

Finally,  knowing  the  desired  thrust  range  and  the  given  pressure  range,  it 
is  possible  to  use  equation  (18)  to  plot  the  thrust  versus  chamber  pressure  for 
various  throat  areas,  as  seen  in  Figure  6.  This  plot  assumes  that  the  value  of  k 
varies  linearly  from  1.432  to  1.783  as  the  pressure  increases  from  0.3  to  10 
MPa.  From  Figure  6,  we  see  that  a  throat  area  of  1  x  1 0*7  m2  will  produce  thrust 


from  0.08  to  1.58  N  while  a  throat  area  cf  1.25  x  10*7  m2  will  produce  a  thrust  of 
0.10  to  1.99  N,  thus  the  thrust  range  for  either  design  is  about  the  same.  These 
particular  values  were  chosen  because  they  come  closest  to  meeting  the  0  to  2 
N  design  requirement,  although  the  throat  area  will  most  likely  be  set  by  the 
maximum  required  thrust,  not  the  minimum  required  thrust. 


Area  Ratio  Vs.  Pressure  Ratio 

for  Nitrogen  (T1  =  -20  deg  C) 


Determining  the  length  associated  with  this  nozzle  requires  the  use  of 
equation  (17).  For  the  1.25  x  10'7  m2  throat  area  and  a  divergence  angle  of 
20°,  the  length  is  3.08  cm,  while  a  divergence  angle  of  5°  yields  a  length  of 
12.81  cm.  For  such  a  disparity  in  length,  one  needs  to  decide  if  the  thrust  losses 
associated  with  a  20°  divergence  angle  is  worth  the  reduction  in  length.  From 
Table  1  we  see  a  thrust  loss  of  3%  for  the  20°  angle  and  1%  for  the  5°  angle, 
both  within  the  error  range  for  the  ideal  nozzle  assumptions.  However,  the 
fragility  of  the  nozzle  and  the  possibility  of  damage  are  probably  much  higher 
for  a  21  cm  nozzle,  thus  a  higher  divergence  angle  is  recommended. 

In  conclusion.  Table  2  summarizes  the  microthruster  design. 


Thrust  Vs.  Chamber  Pressure 


for  Various  Throat  Areas 


Propellant 

Gaseous  Nitrogen 

Specific  Impulse 

72  sec 

Expansion  Ratio 

64:1 

Throat  Area 

1.25  x  10-7  m2 

Chamber  Pressure 

0.3  to  10  MPa 

Chamber  Temperature 

-200  c 

Length 

3.08  cm 

Diuergence  Angle 

20o 

Thrust 

0.1  to  2  N 

Mass  Flow  Rate 

3.4  to  0.1  g/sec 

Table  2:  Microthruster  Design  Parameters 
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Classification:  Repair  (power  distribution)  of  a  partially  functional  asset. 
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Definition  of  Mission  Categories  and  Types 


GTO  Insertion:  perform  GTO  insertion  for  stranded  GEO  satellite. 
Might  also  have  to  perform  subsequent  apogee  (circularization 
and  inclination  change)  maneuver. 
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Replacement:  to  physically  replace  the  dysfunction  component 
with  a  new  unit. 

Supplement:  to  bypass  the  dysfunctional  component  by  supply¬ 
ing  a  new  unit  to  provide  the  function  without  physical  replacement. 
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Mission  Type  Characterization 
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-  Nearly  all  communications  satellites 


Mission  Typs  Chsrsctorizstion 
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Note  that  all  missions  assume  an  automated  rendezvous 
and  docking  capability. 


GEO  Space cra-f-t  Mass  (kg) 


Preliminary  Mission  Analysis 
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Preliminary  Mission  Analysis 
Retrieval 


kit  io  lower  satSiites  to  a  Shuttle-accessible  orbit  for 

pick-up  as  a  targets  of  opportunity. 
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-  Fortunately,  the  propulsive  requirements  are  small  at 
this  distance-even  to  move  a  satellite  out  1,000  km. 

-  Similar  to  a  refueling  mission  except  that  the  kit  would 
use  its  residual  propellant. 
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-  Usable  payloads  to  MEO  will  be  on  the  order  of  100  kg. 

-  Usable  payloads  to  LEO  can  be  nearly  1,000  kg. 
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about  a  20%  performance  increase. 
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